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I. INTRODUCTION 

Radioisotope thermoelectric generators (RTG's) used on unmanned space- 
craft for auxiliary onboard electrical power are commonly fueled prior to launch . 

Unless special provisions are made to modify the environment on the launch pad, 
this means that prior to and during launch there is a period of time that the RTG's 
operate in an air environment. The effects of such air operation on RTG perform- 
ance are basically of two types. First, the air environment temporarily modifies 
RTG performance from its design vacuum values by changing the external and in 
an unsealed system also the internal heat transfer characteristics of the RTG. 

Subsequent vacuum operation is unaffected by these temporary changes in RTG 
operating environment and the RTG will exhibit its original design performance in 
an eventual space environment. Second, operating in an air environment may 
permanently modify the physical characteristics of certain components of the RTG 
such that design performance values are not recovered in a subsequent vacuum 
operation. Mechanical strength characteristics of certain RTG components may 
also be affected. Permanent changes in RTG operating behavior as a result of air 
exposure are usually associated with irreversible chemical and/or physical processes 
caused by reactions and interactions with oxygen. 

Because oxygen poses a serious problem to RTG's that use lead and tellurium 
compounds and alloys , it is accepted practice with devices of this type to hermetically 
seal them. Prior to final sealing, such devices are frequently also back-filled with 
an inert gas in order to reduce the volatility of the thermoelectric material. The result 
of this procedure is that neither permanent nor temporary effects usually occur internal 
to telluride RTG's as a result of air operation. In air, however, the external heat 
transfer characteristics of the RTG are modified by an enhanced heat rejection rate 
from the outer surface of the device because radiation heat transfer is supplemented 
by convection. In air, the RTG consequently operates at hot and cold side temperatures, 
and thus also a performance level, that are lower than the vacuum design values . In 
order for this effect to be only temporary, it is necessary to use materials for the RTG 
casing and radiator, including radiator emittance coating, that are insensitive to oxygen. 
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The air operation of silicon -germanium RTG's presents a situation somewhat 
different from that encountered with telluride generators . Unlike lead and tellurium 
compounds and alloys, silicon-germanium alloys, are not adversely affected by 
operation in air. This is also true of the most common silicon-germanium thermo- 
couple configuration, the so-called Air-Vac thermocouple, that uses components 
which in their respective general ranges of operating temperatures are not adversely 
affected by air, i.e. by oxygen. Inasmuch as volatility is not a problem either with 
silicon-germanium Air-Vac thermocouples, it may appear that silicon-germanium RTG's 
designed for operation in space vacuum require neither hermetic sealing nor any other 
modifications to also permit pre-launch and launch operation in air. Whether this is 
really true or not depends on the sensitivity to air of silicon-germanium RTG com- 
ponents other than the Air-Vac thermocouples . Because relatively little effort has 
been expended to date on the design and development of silicon-germanium RTG's, 
the effects of air operation on RTG performance have never been fully investigated for 
silicon-germanium systems. Nevertheless, because of their proven relative stability 
of performance and reliability, it is namely silicon -germanium Air Vac thermocouples 
that probably will be used in second -generation RTG's intended for use in future 
extended-life space missions. It is the purpose of the present memorandum therefore 
to conduct a preliminary investigation of the effects of air operation on silicon- 
germanium RTG performance. Initially the generator internal components are considered 
exposed to the air environment. After evaluating generator performance under these 
conditions, inert gases are introduced into a hermetically sealed generator and the 
generator performance is re-examined. The gases included in this study are: Argon, 
Helium, Krypton, and Xenon. In terms of the temporary and permanent effects discussed 
above, the primary intent of the present study is to consider the former of these in 
some detail and to point out in a general manner some more important of the possible 
latter-type effects , Because in some applications power is required from the RTG 
even during the launch phase, another and related purpose of the present study is to 
determine the amount of such power available from typical silicon-germanium RTG's in 
air operation and with various internal gases . 
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II 


HEAT TRANSFER MODEL 


The present description of the theoretical aspects of RTG performance in 
both air and vacuum will primarily consider heat transfer effects . The theoretical 
model for the thermoelectric conversion process is independent of generator environ- 
ment and has been throughly discussed previously (see for example References 1 & 2) . 

In vacuum, heat transfer is due only to radiation and conduction; whereas, in air, 
convective heat transfer must also be included . 

This section of the report assumes an idealized generator configuration for 
purposes of a mathematical model. The generator is considered to have radial 
symmetry about an enclosed heat source. The ends of the generator are assumed 
impervious to heat flow (this is equivalent to assuming constant end losses in actual 
RTG calculations) and the sides consist of thermoelements separated by thermal 
insulation. The inner and outer surfaces of the side walls are assumed to be such 
that only radial heat flow occurs within the walls. The hot and cold junction temp- 
eratures of the thermoelements depend on the temperatures of the inner and outer sur- 
faces of the walls; the emissivities of the walls are assumed to be uniform. 

HEAT TRANSFER TO THE RADIATOR 

The total heat transported though the walls of the generator generally consists 
of heat passing through the thermoelements and shunt heat passing through the thermal 
insulation that surrounds the thermoelements. The total heat transferred through the 
walls of the generator is independent of generator environment except for possible modif- 
cation of the thermal conductivity of insulation. In general, the insulation thermal 
conductivity, kj , may be separated into two components; 

k=k + k (2-1) 

I c r 

Where k^ and k^ represent the conductive and radiative components respectively. 
Whereas the conductive term dominates in fibrous insulations, it is the radiative term 
that is usually more important in foil insulations. In general, the two conductivity 
terms have a different power dependence on temperature. For the case of foil insulation 

,,4 * 

it has been found that the two components of thermal conductivity m‘ay be approximated 
by (Ref. 3); 
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k = aT , 
c 

(2-2) 

■k = bT 3 . 

(2-3) 


Where a and b are constants and T is the temperature of the foil insulation. For 
any insulation, however, a power law approximation may be made for the thermal 
conductivity over limited ranges of temperatures as follows: 

k = c T y (2-4) 


As shown in appendix A, an effective (average) thermal conductivity, <k^>, here 
applied to the insulation, may be defined as: 


q I = <k I > 


AT 

l 


= k T (T) 


dT 

dx 


(2-5) 


Where q is the heat flux through the insulation, T is temperature, and X and x are 
length parameters. Integration of Eq. (2-5) yields: 


<k I > 


1 

A T 



(2-6) 


Using average thermoelectric properties (see Appendix A) , the heat balance at the 
hot junction of the thermoelements may be written as: 


Q. =<k >A + <k > A + < k > A — +NI ST 

w in I e X nnX ppX SH 


(2-7) 


2 ' 

- 1/2 NI C R. , 

O 1 
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Where Q. is the total heat input, A is the total effective area of thermal insu- 
m e 

lation, AT is the temperature difference between the hot and cold junction tem- 
peratures, and respectively, JL is the thermoelement length, N is the 
number of thermocouples, I is the series current per thermocouple, R, is the 

U 1 

internal electrical resistance of a couple and A and A are the total cross- 

n P 

sectional areas of the n- and p-type thermoelements respectively. The thermo- 
physical properties of the Seebeck coefficient and thermal conductivity are de- 
noted by S and k respectively. The subscripts n, p and I pertain to n-and p-type 

thermoelements and the thermal insulation. Note that S is defined by S=S -S . 

P n 

Equation ( 2 - 7 ) is one of two basic equations required in the establish- 
ment of thermoelement hot and cold junction temperatures, and thus of thermo- 
electric generator performance . The other such equation pertains to heat rejection 
from the generator and is discussed below for vacuum as well as air operation. 



It is assumed that heat rejection from the radiator is strictly radiative 


in vacuum. The total heat rejected is given by: 


Q , = Q. (I-77) ( 2 - 8 ) 

rej in ' 

where 77 is the conversion efficiency of the thermoelectric generator. Efficiency 
enters the above equation because 77 Q of heat is converted to electrical power 
output. Using the radiation equation, Eq. ( 2 - 8 ) may be rewritten as: 

Q in ( 1 '* )= 0 Vf\ (2 - 9 > 

Where a is the Stefan-Boltzmann constant, is the radiator emissivity, 77 is the 
fin effectiveness, F is the view-factor, A^ is the radiator area, T^ is radiator base 
temperature and T b is the effective external heat sink temperature. 
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After allowances are made for a temperature drop due to conducted 
heat between the radiator base and thermoelement cold junctions, Eqs . (2-7) and 
(2-9) allow the determination of thermoelement hot and cold junction temperatures in 
vacuum operation. Because the efficiency which enters Eq. (2-9), depends upon 
both and T^, it is necessary to solve Eqs . (2-7) and (2-9) iteratively. The 
efficiency is initially given a small but arbitrary value to allow a first approxima- 
tion estimate of T and T,., . The generator performance is then calculated and the 
new efficiency value substituted into equation (2-9) . The iteration is continued 
until the junction temperatures converge to within any desired degree of accuracy. 
HEAT RETECTION FROM THE RADIATOR IN AIR OPERATION 

In addition to radiative heat rejection from the radiator, convective 
heat transfer also takes place in the air operation of a thermoelectric generator. 

For present purposes it is assumed that radiative and convective heat transfer 
terms are additive. Equation (2-9) will accordingly be rewritten as: 

Q in (l-n) - OVUM’S - T ld +M c(V T em 1 / = Q r +Q c ' (2 - 9 10 > 

Where K is the convective heat transfer coefficient, A is the total radiator area 

c 

that contributes to convective heat transfer, 0 is an exponent, and Q and Q 

I* c 

represent radiative and convective heat transfer respectively. The convective heat 
transfer coefficient, K,and the exponent, , depend upon the particular mode of 
convective heat transfer applicable to a given situation. In the case of natural 
convection from vertical cylinders, the numerical values of K and depend upon 
whether turbulent or laminar heat flow exist in the boundry layer between the 
radiator surface and the adjacent environment. It can be shown (Ref. 4) that a 
variable, X, establishes the boundary layer flow conditions in accordance with 
type of flow as follows: 

4 9 

Laminar range 10 <X< 10 

9 12 

10 < X<10 


Turbulent range 



The quantity X may be defined as : 


X = 6.37xI0-k 1 * 3 |T R - T amb 


( 2 - 11 ) 


Where L is the vertical height of generator and 2 is the product of the Prandtl and 
Grashof numbers . 

Figure 1 has been reproduced from Reference 4 to permit calculation of the variable 
Z The film temperature in Figure 1 is defined as the mean temperature between the 
ambient and radiator temperatures . 

Depending on the type of flow, the convective heat transfer from the 
radiator may be written as: 


and 


Q =44.8x10 
c 


-5 


T -1/4' 
xL x 




x[T -T , , 
R a mb] 


1.25 


(2-12) 

(Laminar flow) 


Q = 29.4 x 10~ 5 x A x(T -T A 
c R l R ambj 

Substituting Eqs . (2-12) and (2-13) into Eq. (2-10) yields 


1.33 


(2-13) 

(Turbulent flow) 


Q. (I- 77 ) = ae A fT 4 - T 4 , ]+ 44.8 xlO 5 x L 1//4 (2-14) 

in 11 r Rl R ambj 

1.25 (Laminar flow) 


x A T _ -T 
R C amb 


and 


Q in (1 ^> = -Cbl + 29 ' 4 * 10 ' 5 X A 


R 


(2-15) 


1 « 

(T -T A (Turbulent flow) 

R ambj 

where the fin effectiveness and view factor are taken to be unity and the convective 

and radiative area are considered equal (A = A ) . 

C Ja. 
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Equations ( 2 - 7 ), (2-14), and (2-15) are the equations necessary to 

determine the thermoelement hot and cold junction temperatures in air operation . 

As before, the dependence of Eqs . (2-14) and (2-15) on efficiency again requires 

an iterative calculation . The non-integral power dependence of temperatures in 

Eqs. (2-14) and (2-15) introduces an additional complication. The Newton- 

Ralphson approximation technique may be used to determine T . A discussion of 

K 

this technique may be found in most books on Advanced Calculus . 
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III. COMPUTATIONAL METHODS 

In order to calculate the performance of a thermoelectric generator in 
air as well as in vacuum, the dimensional parameters of the overall generator 
and generator components are required as inputs to the model discussed in the 
preceding section. Thermophysical properties of the thermoelements, thermal 
insulation and heat rejection system are determined from previously reported 
data. With this information and the known heat input, the performance of any 
thermoelectric generator may be established. The only uncertainties which may 
exist in the performance calculations are those associated with the thermophysical 
properties of the various materials used and the heat transfer model employed. 

In the following section, in which actual performance data are compared to per- 
formance calculated with the model described in this memorandum, it will be 
seen that close agreement is obtained between theory and measurement. In cases 
where performance data exist and precise thermophysical property data are lacking, 
the described computational methods may be used to back calculate the unknown 
property data . 

HEAT BALANCING 

As a starting point in the calculation of thermoelectric generator per- 
formance, the heat must be balanced for the system. Since the performance of 
the power system is governed by heat emanating from the heat source and the re- 
sultant temperature gradients, a reasonable estimate of various losses must be 
made. The heat source is considered to be completely surrounded by either 
thermoelectric elements, thermal insulation or heat source support members. Al- 
though the exact configuration of the generator is not of particular importance, 
careful consideration of effective areas must be made , Figure 2 depicts a typical 
cylindrical thermoelectric generator and illustrates some important considerations 
in the balancing of heat . 

Heat flowing from the heat source may be decomposed into radial and 
axial components . The axial flow is composed of heat conducted through the heat 
source support structure and that lost through thermal insulation on the ends . These 
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heat losses serve no useful function in the energy conversion process; quite the 
contrary, they detract from the heat available to the thermoelements . The radial 
heat is similarly divided into two components . The first component comprises the 
heat that passes through the thermoelements, while the other component consists 
of heat lost through the thermal insulation . Both of these radial heat flow com- 
ponents are considered to consist of simple conduction processes that may be 
described by means of effective values of thermal conductivity and cross-sectional 
area. The thermal conductivity values are integrated averages over the appropriate 
temperature ranges being considered . 

To determine the radial heat flux passing through the thermal insulation, 
an effective insulation area must be used. Since this area is dependent upon 
radius, the radial heat flux density may be written as: 


= _Q_ = 

2nrL 


k 

K dr 


(3-1) 


Where Q is the total radial heat flowing through the insulation, r is the radius 
and L is the axial length of insulation. 

Rearrangement and integration of Eq. (3-1) yields: 


r=b 

q[ 


2ttL <k> AT, 


(3-2) 


Or: 


Q = ~~- < k > AT. (3-3) 

Jinb 

a 

Where < k> is the effective thermal conductivity of the insulation, b is the outer 
radius, and a is the inner radius of insulation. 

A comparison of Eq. (3-3) with the conductive heat transfer equation, allows the 

definition of the effective cross-sectional area, A ,, , or radius, r , as: 

eff eff 
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(3-4) 


A eff = 2TTL fe = 2T,r eff L • 

' a 

For inner and outer radii typically encountered with thermoelectric generators, the 
effective radius is very nearly b+a , i . e . the linear average of the inner. and outer 
radii . 

The conductive temperature drop between the radiator base and thermo- 
element cold junctions may be calculated by use of dimensions and properties of 
appropriate components and the total heat being transported. A reasonably good 
estimate of this temperature drop jAT^j may usually be made . Temperature drops 
of the order of 20 to 40° C are typically calculated. 

Not all of the heat rejected from a thermoelectric generator is transported 
by radiation and/or convection. Since the generator must be mechanically supported, 
a fraction, cct of the rejected heat is carried off by conduction through structural 
members . Usually this is only a small fraction of the totalrejected heat, but for 
accuracy has to be taken into account. Similarly, for accuracy, it is also neces- 
sary to account for radiator emissivity-, e , and fin effectiveness, 71 , and the 

r 'f 

radiation view factor F . 

The heat rejected by the radiator, Eq. (2-10), may be rewritten as 

Q rej - (l-J?-a)Q in - VW( T R "O + M c ( T R - T a (3 - 5> 

Where all of the symbols have been previously defined. 

The calculation of the effective convective and radiative areas of a 

thermoelectric generator, A^ and A R respectively, is best illustrated by reference 

to Figure 2 . The effective convective area, A , may in the first approximation 

c 

be assumed to be given by the total external surface area of the generator. The 
effective radiative area,A^F, may to the same approximation be assumed equiva- 
lent to the total equivalent area of the generator indicated by the dashed lines in 
Figure 2. This represents the minimum surface area required to enclose the RTG. 

The quantities a. and in Eq. (3-5) have also to be determined or approximated 

prior to the use of the equation to calculate the radiator base temperature T . 

K 
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As already mentioned, Eq. (3-5) is solved by the Newton-Ralphson approximation. 

Ignoring Thomson heat, the heat balance at the hot junctions of the 
thermoelements may be written as 


Q =A.Q. = [<k T >A + N<k > a +N<k >a (3-6) 

radial in lie n n p pJ « 


+ NI g ST H - l/2NIg 6 


<P n > <o p > , 

“sl + a 
n p 


Where a and a are the cross-sectional areas of individual n- and p-type thermo- 
n p 

elements, X is fraction of total input heat that traverses the generator in the radial 
direction, 5 accounts for extraneous circuit resistance and p^ and p^ are the 
electrical resistivities of the n-and p-type thermoelements . The definition of a 
total radial thermal conductivity, K , allows Eq. (3-6) to be rewritten as 


«in = V T+ (^V 1 + 2 m)+T c], 


(3-7) 


where 


v _ N 
x j e 


<k T > -rf- +<k > a + <k >a 
IN n n p p 


m = R jjR g 


r = 


NS 


26 tn + h> 

l a n a p- 


Note that R T is the external load resistance and R is the internal generator resist- 

L g 

ance . Using the radiator temperature determined by Eq. (3-5) and recalling that 


T - T + AT 
C k X C 


(3-8) 


allows solution of Eq. (3-7) for T . 

H 
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Equations (3-5), (3-7), and (3-8) have been used to establish both the 
hot and cold junction temperatures . The non -integral power (3 and the efficiency 
term (l~rj) in Eq. (3-5) require two discrete iteration schemes in order to deter- 
mine the temperatures of the generator. Because generator temperatures are not 
known a priori, it is necessary to start the calculational sequence by using 
arbitrary values of thermoelectric properties . As new temperatures are calculated 
in each iteration, an integration scheme is used' to calculate average thermo- 
electric properties appropriate to these temperatures in each iteration. 
PERFORMANCE CALCUIATIONS 

Once the thermoelement junction temperatures have been established, 
the performance of the generator is calculated in a straight-forward manner . The 
load voltage is given by: 



m 

(1+m) 


N a SAT . 


The power output may be written as: 


(3-9) 


P 


m 

out (l +m) 2 


2 2 
NS AT 

6R. 


The load current is: 


(3-10) 


j =n-- SAT 
X L 6R, 


(3-11) 


Where m is the load ratio (R^/R^) and n is the number of parallel electrical branches 
in the generator. Generator conversion efficiency is given by: 


n = P Jq. (3-12) 

1 out / in 

Because of the iterative nature of the calculational sequence outlined 
above, it is particularly convenient to effect the solution on a high speed computer. 
An example of the corresponding computer program, actually used to obtain the re- 
sults discussed in Sections IV and V of this memorandum is reproduced in Appendix 

B. 
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IV. RCA MODULE EVALUATION 


In order to test the validity of the mathematical model discussed in Sections 
H and III and, in some instances, to calculate thermal conductivities of selected 
thermal insulations , this section of the present memorandum deals with the calcula- 
tion of the performance of RCA converters C-2 , C-3 and C3Z, The calculated 
performance values are compared to the experimental data reported in References 6 and 
7. The thermoelectric material used by all of these converters is a silicon-germanium 
alloy of composition 63 a/o Si - 37 a/o Ge in the conventional "Air-Vac" configuration 

with silicon -molybdenum alloy hot shoes . Both the n-and p-type thermoelements 
2 

have a 0.0889 cm cross-sectional area with an active length of 2.155 cm in all 

three converters.. All converters employ eighteen thermocouples arranged electrically 

in series and have an hexagonal overall geometry . The thermal insulation in each 

2 

converter has an effective area of 171 cm . 

The temperature drop between the cold junction of the thermoelements and the 
radiator base is calculated to be approximately equal to 0.5 times the heat input. 

This figure is based upon approximate thermal resistances of the various components 
of the cold stack and the connection to the radiator surface. The primary difference 
between the three converters in question is the type of thermal insulation and the exact 
radiator configuration of each. 


CONVERTER C-2 

The C-2 converter employs multiple radiation/convection fins for waste heat 

rejection and uses Min-K 2020 thermal insulation. The radiation area of the radiator 

2 

is calculated to be 817 cm , while the effective convective area is estimated to be 
2 

about 1200 cm . The heat rejected from the radiator in vacuum is given by Eq. (3-5) 

with K set equal to zero. Using the RCA experimental data and the radiative area 

indicated above, Eq. (3-5) may be solved to give e 7? , F = 0.76. This value may be 

r f 

expected to remain nearly the same for all the RCA generators considered in this report. 
Exact agreement in the cold side temperature (T _ ) , over the entire input power 

range, .is noted in Fig. 3 for both vacuum and nitrogen environments. (Exact agreement 
is, of course, necessary at the P IN =56 WATT vacuum point since the corresponding 
temperature was used to establish the value of F = 0.76.) In operation in an air 
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(nitrogen) environment, the convective heat flow is laminar and the gemetry 
dependent convective heat transfer coefficient, K, is calculated to be 
2.38 x 10 ^ watt/cni - °K. The fraction, 1 - X , of heat lost from converter 
ends was assumed to be 0.05. 

An effective value of thermal conductivity for the thermal insulation 
must be used in the calculation of thermoelement hot junction temperatures. 

The thermal conductivity of Min-K 2 020 reported in Ref. 5 was used as a starting 
point. Inasmuch as the authors of the referenced report, express doubt about 
the validity of the reported values, these values of thermal conductivity were 
only used in the initial calculations and later adjusted to yield agreement between 
calculated and experimentally determined performance of the C-2 converter. As 
a result, the effective thermal conductivity of Min-K 2020 was found to be 
slightly greater in vacuum and slightly less in air than the values reported in 
Ref. 5 . As seen in Fig. 3, exact agreement between the calculated and experi- 
mental hot junction temperatures was not achieved at the higher power input 
levels . Agreement between calculated and experimentally determined performance 
for the converter is however excellent (Fig. 4) . This agreement is not surprising 
since it was the experimental power output of the converter that was used to back- 
calculate the values of effective thermal conductivity for Min-K 2020 . Similar 
discrepancies in the thermoelement hot-junction temperatures were found to exist 
in the performances calculated for all of the RCA converters; this in spite of the 
fact that in some cases values of thermal conductivity back calculated by RCA 
were used. The small discrepancy between calculated and measured hot junction 
temperatures may possibly be partly attributed to the fact that the calculated tem- 
peratures refer to the actual hot-junctions whereas the measured temperatures 
most likely are dependent upon the location of the temperature sensors on the hot 
shoes . A small discrepancy in the thermophysical properties of the thermoelements 
may also result in some difference between calculated and experimental tempera- 
ture values . 
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CONVERTER 03 

The C-3 converter employs six radiation fins for waste heat rejection and 

2 2 

uses multi-foil thermal insulation. Effective areas of 719 cm and 952 cm were 
calculated respectively for the effective radiative and convective radiator areas. 

The product e r) . F was again taken to be 0.76 and the convective heat transfer 

r f 4 2 o 

coefficient was again found to be 2.38 x 10 watt/ cm - K. Converter C-3 is 

identical to C-2 except for the radiation/convection fins and the thermal insula- 
tion of the converter. The effective thermal conductivity of the multi-foil insula- 
tion in both nitrogen and vacuum was obtained from Ref. 6. Calculated thermo- 
element temperatures are shown in Fig. 5. Experimental data points measured at 
RCA are included for comparison. Other than the slight discrepancy in hot-junction 
temperatures which has already been mentioned, good agreement exist between 
calculated and measured data. Fig. 6 shows plots of converter power output and 
efficiency as functions of heat input. Good agreement between calculated and 
experimental data is again observed. 

CONVERTER C-3Z 

The C-3Z converter is identical to the C-3 except for the thermal insulation. 
Whereas converter C-3 employs multi-foil insulation, converter C-3Z uses Zircar - 
(ZrOJ instead. The absence of reported information on the effective thermal conduc- 
tivity of Zircar in nitrogen necessitated a back-calculation from experimental 
performance data . Figs . 7 and 8 show plots of calculated thermoelement hot and 
cold junction temperatures and power output and efficiency for the C-3Z converter 
as functions of power input to the converter. As before, experimental data reported 
by RCA (Ref .7) are shown for comparison. The agreement between calculated and 
experimental data are once again on the whole quite good. 

For convenience the back-calculated effective thermal conductivity values 
for the multi-foil, Zircar and Min-K 2020 thermal insulations in both nitrogen (air) and 
vacuum are plotted in Fig. 9. As already discussed, some of these back-calculated 
thermal conductivity data were obtained from RCA and some were an outcome of the 
present work. 
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V. MULTI-HUNDRED WATT GENERATOR EVALUATION 

In evaluating the performance of the MHW-RTG, several distinct cases will 
be considered. Inasmuch as some question still exists on the type of thermal 
insulation to be used in the MHW generator, the present calculations have been 
performed for all three types of insulation considered in Section IV: the multi-foil, 
Zircar and Min-K 2020. The first case (Case I) is concerned with a generator 
geometry designed to optimize a RTG employing foil insulation which requires a 
thermoelement length of 0.8 inch. The generator is initially considered "open", 
thus, allowing direct exposure of the internal components when operated in air. 

In the second case (Case II), the generator has an approximate 18% increase in 
the thermoelement £/A ratio to permit RTG performance optimization with 
fibrous insulation. Again the generator is initially considered "open". The 
generator configuration and dimensions used in the present analysis are those 
reported by GE at the MHW design review meeting, 10 December 1969, Washington 
D.C. The more important design parameters are tabulated for reference in Table 1 . 

By simply modifying the basic MHW design, it is possible to hermetically seal 
the RTG such that an inert gas may be contained within the generator during the 
air-operation phase. As will be seen, this eliminates many of the problems associated 
with generator operation in air. Although the generator exterior operates in air, the 
temperatures are low enough to avoid most problems associated with chemical reaction, 
etc. By selecting a low thermal conductivity gas, such as Xenon, the generator will 
deliver a significant amount of power in this "pre-launch" condition. After reaching 
the eventual vacuum environment of space , the generator could be vented to allow 
its performance level to be re-established at that corresponding to vacuum operation. 

The use of inert gases in this application is presented by first considering the 
thermal conductivities of Argon, Helium, Krypton, and Xenon. The effect of these 
gases on the conductivity of insulation is then considered and finally the 
MHW-RTG performance is re-examined and a summary is provided (Tables II and III). 
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Case I Thermoelement Length. = 0.8 inches 

The generator is initially considered "open" and makes use of a silicon- 
germanium alloy with composition 80 a/o Si - 20 a/o Ge in the standard Air-Vac 
thermocouple configuration. The generator is a hexagon with conical end-caps. 

A schematic of the generator is shown in .Fig. 10. The view factor, F, of the 
radiator has been assumed to be unity and the fin effectiveness has been set at 
0.95 (Ref. 8) in the present analysis. It has been further assumed that 2.5 percent 
of the total heat input to the generator is lost through each end-cap; the remainder 
of the heat, less the power output, is rejected by radiation and/or convection. 

The radiator area, the same in this case for both radiative and convective heat 
transfer, includes the total lateral surface area of the generator and of one end-cap. 
The radiation area of the second end-cap has been discounted because in its 
intended application the generator will be mounted axially in tandem with another 
generator. A temperature differential of about 20°C has been calculated between 
thermoelement cold junctions and the radiator base by accounting for the thermal 
resistances of all cold stack members and the heat passed through each stack. 

Using RCA supplied thermoelectric property data for the 80 a/o Si - 20 a/o Ge 
alloy, initial calculations of the foil insulated MHW generator performance in 
vacuum were conducted to test the applicability of the model. The results of this 
initial calculation are illustrated in Fig .11. Calculated data reported by GE are 
shown as data points in the figure. Agreement between the performance calculated 
with the model described in this memo (using RCA property data) and the calculated 
performance values reported by GE is good. Although not shown in Fig. 11, similarly 
good agreement was obtained for the thermocouple cold junction temperatures. 

In view of certain questions regarding the thermoelectric property data 

reported by RCA for the 80 a/o Si - 20 Ge alloy, the calculation was 

repeated with the "1500 hour" property data of the alloy given in Memorandum #6 

(Ref. 9). The results of that calculation are shown by the dashed curves in Fig. 11. 

It is noted that the agreement between GE data and the present data are no longer as 
good as before. The difference in the two sets of data is however, primarily due to 
the material properties and thus has little to do with the model being used. In fact. 
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the validity of the model has again been established by the close agreement of 
the GE data with those presently calculated when using identical thermoelectric 
properties. For purposes of consistency, however, all subsequent MHW 
generator performance calculations used the thermoelectric property data reported 
for the 80 a/o Si - 20 a/o Ge alloy in Ref. 9. 

The hot and cold junction temperatures and performance of the multi -foil 
insulated MHW generator (thermoelement length = 0.8 inch) in air and vacuum 
operation are shown as functions of load current in Figs. 12 and 13. Generator 
performance in -vacuum has assumed ambient heat sink temperatures of 0 and 311°K. 
Only the latter heat sink temperature has been assumed for the air operation case. 
Although generator operating temperatures are dependent on the ambient heat sink 
temperature (see Fig. 12), the performance of the generator is practically indepen- 
dent of it. Thermoelectric properties and temperature-differential across the 
thermoelements are minimally affected by small changes in generator radiator 
temperatures . The power output and load voltage as functions of load current 
(see Fig. 13) are the same for the 0 and 311°K ambient heat sink temperatures. In 
addition to radiation, the air operation case has assumed heat rejection by natural 
convection; for the conditions of the MHW generator, the convective heat flow has 
been calculated to be laminar. 

The hot and cold junction temperatures and performance of the Zircar-insulated 
MHW generator (thermoelement length = 0.8 inch) are shown as functions of load 
current in Figs. 14 and 15 for the same ambient heat sink conditions as assumed in 
connection with the foil-insulated generator. The corresponding data for the generator 
using Min-K 2020 thermal insulation are plotted in Figs. 16 and 17. 

Case II Thermoelement Length =1.2 inches 

The results of the calculations performed with an increased thermoelement length 
of 1.2 inches ( l/h increased - 18%) are shown in Figs. 18 through 23 for the three 
different insulation materials. As indicated in Fig. 18, excessively high hot -junction 
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temperatures ( ~ 1200°C) result when the increased thermoelement length is used 
in conjunction with foil insulation. These high temperatures discount using foil - 
insulation in this case . 

On the other hand, the use of Zircar insulation with the increased length has 
resulted in performance (see Figs. 20 and 21) comparable to that calculated for 
foil-insulation with the smaller i/h ratio. With the Zircar-insulatibtr, maximum 
power is observed at a hot -junction temperature of 1125°C. In view of some recent 
experiments, this temperature is probably excessive for Zircar and some 
insulation decomposition into metallic components would result. 

Figures 22 and 23 show the calculated performance of the increased thermo- 
element length generator using Min-K 2020 insulation. Improved performance is 
noted; but again, high temperature problems would most probably be encountered in 
vacuum operation because of excessive insulation temperatures. 

The thermal conductivities of the insulations used in the calculation of the 
MHW performance are those plotted in Fig. 9. It should be noted that the air- 
operation calculations have thus made use of effective insulation thermal conductivity 
values in nitrogen. Although insulation thermal conductivity is probably slightly 
different in air than it is in nitrogen, the difference is not expected to be great. It 
should also be noted that the thermal conductivity of Min-K 2020 has to be extrapolated 
over a large temperature range (up to 1100°C) before it can be used in the vacuum 
calculations . This may introduce some error in the calculations if the conductivity 
deviates greatly from the projected behavior. The MHW generator operating temper- 
atures and performances shown in Figs. 11 to 23 for air operation are however, 
believed to be reasonably accurate. 

It should be noted that the performance of the "open" geometry MHW generator 
in going from vacuum to air operation is significantly reduced for all three thermal 
insulations considered. The biggest performance reduction occurs in the case of the 
multi-foil insulated generator, with the smallest reduction for the case of the 
generator that uses Min-K 2020 insulation. Vacuum performance, however, is greatest 
for the case of the multi -foil insulation and least for the Min-K 2020. 
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General Discussion of "Open" Generator Configuration 

Of the two cases considered, the generator using the shorter thermoelement 
lengths (Case I) appears to give the most desirable overall performance. The high 
hot-junction temperatures encountered in -'Case II, would very likely cause insula- 
tion problems when the generator is operated in vacuum. Consequentially, the 
following discussion will be primarily directed towards the generator employing a 
thermoelement length of 0.8 inch. Whereas Figs. 12 to 17 show the relative 
temperatures and performances of the Case I MHW generator in air - and vacuum - 
operation for three different thermal insulations, the data in these Figures of course 
give no direct indication of potential problems that may arise as a result of the 
operating temperatures/environments in terms of damage to the generator. Considering 
first the case of the multi-foil insulated generator, it is noted that under open 
circuit operating conditions in vacuum the hot junction temperatures of the thermo- 
elements approach the solidus temperature of the 80 a/o Si - 20 a/o Ge alloy 
(about 1300°C) . The fuel capsule temperature will in this case be of course 
correspondingly higher. It should be obvious that such a condition is undersirable 
from the standpoint of possible permanent damage due to a variety of effects 
associated with sublimation, chemical reactions , melting, sintering, etc. Every 
precaution should therefore be exercised to prevent the MHW generator from ever 
operating in an open circuit operating mode. In air operation the foil insulated MHW 
generator will exhibit thermoelement hot junction temperatures of the order of 800 to 
900°C. As regards silicon-germanium Air-Vac thermocouples, air operation at such 
temperatures is permissible without any deleterious effects . The cold side bonds 
to tungsten electrodes , as well as all other cold side components , should be able 
to withstand air operation for reasonable lengths of time at the indicated cold junction 
operating temperatures. Possible problems that may be encountered with the foil- 
insulated MHW generator are thus generally related to generator components other 
than those associated with the Air-Vac thermocouples. Probably the two main areas 
in which problems may occur are the multi-foil insulation itself and the graphite 
re-entry shell that encapsulates the fuel capsule in a RTG. At the indicated operating 
temperatures, refractory metals will severely oxidize in air and thus should not be 
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used in an unsealed MHW generator. Possible materials that may be used for 
the foil insulation if it is to operate in an air environment at elevated temperatures 
are nickel, platinum and various other nobel metals. Unfortunately, however, 
may of these metals are not compatible, with silicon-germanium alloys at elevated 
temperatures because of relatively low temperature eutectic formations . Thus , for 
most non-refractory metal foil insulations to be. used in silicon -germanium generators 
that operate at elevated temperatures, it is necessary to replace a portion of the 
■inner foil insulation with fibrous-type insulations that are compatible with silicon- 
germanium alloys and that are designed for high temperature use. Although refractory 
metal foil insulations probably are compatible with silicon-germanium alloys at even 
the highest temperatures at which such alloys can be used, their sensitivity to 
oxidation precludes their use from unsealed silicon-germanium RTG's that will undergo 
air operation. In conclusion then, for refractory metal foil insulations to be used in 
high temperature silicon-germanium thermoelectric generators it is necessary that 
the inside of the system not be exposed to air during generator operation. If it is 
desired to operate an unsealed generator in air, it will be necessary to use non- 
refractory metal foil insulation or fibrous insulation. In the former case it will be 
necessary to use hybrid fibrous-foil insulation system to eliminate possible 
problems of interaction of the foil insulation with silicon-germanium alloys . As regards 
the use of a graphite re-entry shell to encapsulate the fuel capsule of a RTG, 
oxidation problems very definitely would occur at the indicated hot side operating 
temperatures if the system were to. operate in air. The system would either have to be 
sealed or the graphite shell protected with an oxidation resistant coating in air 
operation. 

As regards the fibrous insulations, Zicar (Zr0 2 ) and Min-K 2020 con- 
sidered in the present analysis of MHW generator performance, it is noted that 
at the maximum power output operating point they both operate at temperatures com- 
patible with their high temperature capabilities in both air and vacuum for the Case I 
thermoelement length of 0.8 inch. The case of the Min-K 2020 insulation in 
vacuum, however, is somewhat marginal in that its highest recommended 
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operating temperature is of the order of 1000°C. Under open circuit conditions 
with the MHW generator hot side vacuum operating temperatures in the 1100 to 
1200°C range, it is believed that severe degradation problems will probably be 
encountered with both types of insulation. This condition will be even worse if 
the fibrous -insulated MHW generator is re-designed to operate at the maximum 
power point with a hot junction temperature of 1100°C (as in Case II) . The 
points made in regards to possible oxidation problems of the graphite re-entry 
shell in air operation are also valid in the case of the fibrous insulated MHW 
generator. As before, no special problems are to be anticipated with the 
silicon-germanium Air-Vac thermocouples in either air or vacuum operation. 
Enclosed Generator Employing Inert Gases 

The introduction of a protective gas into the internal portion of the MHW 

generator reduces many of the forementioned problems associated with generator 

* 

operation in an oxidizing environment. The analysis of this mode of operation 
follows that given in Ref. 10. The essential details of the analysis are discussed 
in the remainder of this section of the memorandum. To perform the analysis the 
thermal conductivity of Argon, Helium, Krypton and Xenon was determined, then the 
effect of the gas thermal conductivity upon the converter insulation was evaluated 
and finally the performance of the RTG was calculated. 

Thermal Conductivity of Gases 

The thermal conductivity of several gases of interest are shown in Fig. 24 . 
The solid lines represent experimental data for several of the gases . The dashed 
lines represent the calculated thermal conductivity using the technique given in 
Ref. 11. This technique used the Sutherland formula for computing the viscosity, 

/i , of a gas as: 


5_ ( K MT) 1/2 

^ 16 or 1+S_ 

T 


(5-1) 


Where ^ is collision diameter of gas molecules, K is Boltzmann's constant, M 
is molecular mass , T is temperature, and S is the Sutherland's constant. 
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When the viscosity at one temperature , T , is known the Sutherland formula 

o 

can be rewritten as 


3/2 

a = u > T x ' ,T + S v 0 v 

p o (-) (_o ) (5-2) 

o T + S 

k 

The product, — — , where k is the conductivity and C is the specific heat 

fj, O 

at constant volume, is practically independent of temperature. Thus for gases 
with a constant specific heat the conductivity varies as the viscosity. The 
difference between the calculated values and experimental values for air and 
nitrogen reflectsthe fact that the specific heat does vary somewhat with temp- 
erature, and that nitrogen is diatomic rather than monatomic, and in the case of 
air, it is a mixture rather than a pure gas. However, in the case of other gases 
good agreement between the calculated and experimental data is achieved (Ref. 12 
and 13) . All of the gas conductivity values are considered valid for pressures of 
about 0.1 to 2.0 atm . 

Effect of Gas on Insulation Conductivity 

As has been shown, the conductivity of various insulations of interest 
increase when used in a gas environment as opposed to a vacuum environment. The 
total conductivity of a material (insulator or otherwise) can be expressed as: 

K = k + k + k (5-3) 

c r con 

Where k is the contribution due to solid conduction , k is the contribution due to 
c r 

radiation, and k is the contribution due to convection. If the material is operated 
con 

in a vacuum, then: 

K = k = k + k (5-4) 

v c r 

Where k^ is the effective' conductivity in a vacuum. If a gas is then added, the 
conductivity becomes: 

K = k + Ak (5-5) 

v g 
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Where k^ = conductivity of the gas, and A = empirical constant. The value of 
the empirical constant A reflects the degree to which the conductivity of the gas 
is simple additive, i.e. superposition theory is valid, and the degree to which 
convective heat transfer may be occurring in the gas . Thus the value of A would 
depend upon the density of the insulation and the amount of open space within it. 

Fdr example , in an insulation with a high pore volume , A would approach 1 . 0 if 
no convection is occurring . In a dense material with little or no pore volume , A 
would approach zero. In some cases where significant pore volume is present 
convection can be significant. In these cases the value of A could exceed 1.0. 

Using Eq. (5-5) to evaluate A for various types of insulations based upon 
the curves shown in Fig. 9 and. Ref. 14, it is found that for foil insulation, A 
varies from 1.0 to about 1.4. For Zircar, A varies from 0.7 to 1.0 and for Min-K , 

A varies from 0.4 to 0.75. Thus it appears that a significant amount of convection 
does occur in multifoil insulations whereas in the fibrous types the effect of con- 
vection is less severe. The high A value for foil insulation may also reflect some 
changes in the surface emissivity as the gas is absorbed ( i.e.: constant). 

This would be especially true in an oxidizing gas such as air. It should also be 
kept in;mind that the use of Eq. (5-5) to evaluate the conductivity of insulations 
in gas requires the knowledge of the conductivity of the insulation in vacuum . As 
the conductivity of foil and fibrous type insulations are highly dependent upon edge 
effects in vacuum operation, these uncertainties will be reflected into the uncertainties 
of the conductivity when used in a gas . 

Using average values for A,the effect of the gas conductivity upon the insulation 
conductivity has been plotted in Fig. 25. Notice that the effectiveness of the three 
insulations is reversed when operating them in air as opposed to vacuum. Operation 
in Krypton results in about the same conductivity for all three insulations. The points 
at which the conductivity of the gases are noted in Fig. 25 are based upon an average 
gas temperature typical of that encountered in the MHW generator employing thermo- 
element lengths of 0 . 8 inch (Case I) . If the longer thermoelements were used (1 .2 inch) , 
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the average insulation temperatures would increase and as a result, the curves 
would be shifted upward by a significant amount. This would be particularly 
true in the case of Min-K, since the vacuum thermal conductivity is strongly 
temperature dependent. This increase in conductivity has been included in the 
performance calculations . 

MHW Generator Performance Using Inert Gases 

To obtain generality, the electrical performance of the MHW generator has 
been evaluated as a function of insulation thermal conductivity. For the two 
generators considered. Case I - thermoelement length = 0.8 inch and Case II - 
thermoelement length = 1.2 inch, the performance is shown in Fig. 26. The sink 
temperature was held at 311°K for all cases investigated. However, as has been 
shown, the power output does not significantly depend upon the sink temperature 
and therefore the output power can be shown as a function of insulation conductivity 
alone . The behavior is illustrated in the figure . Shown are the peak power and the 
power at a load voltage of 30 volts. Note that the power at 30 volts coincides with 
the peak power for the 0.8 inch element lengths (Case I) at insulation conductivities 
corresponding generally to operation in vacuum; whereas, the 30 volt power coincides 
with peak power for element lengths of 1.2 inch at conductivity values corresponding 
to operation in air. 

Just as the power output is independent of sink temperature it is also inde- 
pendent of RTG external environment, This is not true for the RTG temperatures . 

The temperatures depend both upon the sink temperature and the presence or absence 
of gases external to the RTG. For all cases evaluated using an internal vacuum 
environment , the external environment was also considered to be a vacuum . For 
all cases using a gaseous internal environment, air was assumed to be the external 
environment. The gas pressures were assumed to be about 1 atm. 

Fig. 27 shows the hot and cold junction RTG temperatures as a function of 
insulation conductivity. The temperatures at open circuit, peak power, and short 
circuit are identified. The data shown in Fig. 27 are for a thermoelement length 
of 0.8 inch (Case I) . Fig. 28 shows the temperatures for the case of a thermoelement 
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length of 1.2 inch (Case II). As would be expected the hot junction temperatures 
for the longer element length are significantly higher as the length to area ratio is 
about 18% higher for that case. The cold junction temperatures are essentially the 
same for both element geometries. 

The hot junction temperature at a load voltage of 30 volts may be either 
higher or lower than the temperature at peak power depending upon the particular 
insulation conductivity and element geometry. As shown in Figs. 29 and 30 the 
30 volt power may lie either to the short circuit or open circuit side of the peak 
power current. In general, for the 0.8 inch element length, the 30 volt power lies 
to the open circuit side of peak power for the high conductivity insulation; for the 
1.2 inch element length, the 30 volt power points lie to the short circuit side of 
peak power. Thus for the 0.8 inch element length, the 30 volt hot junction temper- 
atures may be somewhat greater than the peak power temperatures and for 1.2 inch 
element length the 30 volt hot junction temperatures will probably be less than the 
peak power temperatures . 

A summary of the peak power output, 30 volt power output, hot and cold 
junction temperatures at peak power is given in Tables II and III for the various 
insulations and gases that were considered in this study. In vacuum and low 
conductivity gases, foil gives the best performance and Min-K the poorest. For 
air and other gases of comparable conductivity Min-K gives the best performance 
while the foil gives the poorest. For all cases Zircar gives performance in between 
foil and Min-K. See pp. 19 and 20 for comments regarding RTG operating temperatures. 

Two points worthy of further consideration should be mentioned. First, as was 
pointed out , the vacuum performance of foil and fibrous insulations is highly dependent 
upon end losses and thus the geometry of the converter. Consequently, a detailed • 
analysis of the MHW-RTG geometry with respect to this problem needs to be carried 
out to ascertain the validity of the vacuum values of conductivity used in this study. 
Secondly, the computer program used in this study considered the RTG end losses to 
be fixed. In general, if the RTG is operated in a gas the end losses will increase 
by an amount dependent upon the conductivity of the gas . Thus the power output and 
the hot junction temperatures will be somewhat lower than shown here. The magnitude 
of this difference needs to be examined. 
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SUMMARY 

The present study has been concerned with air/inert gas operation of thermo- 
electric generators designed for vacuum operation. The object of the study has been 
the development of a model to determine gas and vacuum performance of such generators. 
The validity of the model has been verified by its applications to several existing 
generators which have been tested in both vacuum and air (actually nitrogen) environ- 
ments . The model has also been used to back- calculate the thermal conductivities 
of selected thermal insulations that presently are considered for use in high temper- 
ature silicon-germanium RTG's. Finally the model has been applied to the MHW 
generator, with two different thermoelement lengths, to determine its' performance 
under various environmental conditions. The initial characterization was performed 
in vacuum and then the effects of various inert gases upon the insulation conductivity 
were included . The calculations were performed with three different thermal insula- 
tions: multifoil, Zircar and Min-K 2020. Conclusions drawn from the latter evaluation 
may be summarized as follows: 

1. The foil insulated MHW generator employing 0.8 inch thermoelements 
exhibits some 54 percent of its peak-power vacuum performance in air 
operation. Use of 1.2 inch thermoelements in foil- insulated generator 
results in excessive temperatures. The MHW generator insulated with 
Zircar exhibits some 59' percent for Case I and 75 percent for Case II of 
its peak-power vacuum performance in air operation; while, the 
•corresponding percentages, for the Min-K 2020 insulated generator are 
73 percent and 89 percent for the Case I and Case II respectively. 

2. For fixed generator design, foil insulation yields the best peak-power 
vacuum performance for the MHW generator. Zircar insulation gives an 
intermediate performance level; while, Min-K insulation gives the lowest 
performance. 
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3. No permanent damage to the Air-Vac thermocouples used in the 
MHW generator is anticipated as a result of air operation with any 
of the internal gases considered. 

4. If the MHW generator is permitted to operate in an open-circuit 
operating mode in vacuum, excessive hot side operating temperatures 
will probably cause permanent damage to the Air-Vac thermocouples 
and also the thermal insulation and possibly to the fuel capsule. 

5. The use of foil insulations in the MHW generator requires special 
consideration when planning to have an "open" generator operate 
both in air and vacuum . Refractory metal foils are suitable for vacuum 
operation but cannot be used in -an unsealed generator that also 
operates in air. Nickel, platinum and other noble metal foils can 
operate in both air and vacuum but at elevated temperatures are incom- 
patible with silicon-germanium alloys because of relatively low tempera- 
ture eutectic formations . 

6. No incompatibility problems exist with the fibrous insulations 
(Zircar and Min-K 2020) up to hot side operating temperatures of the 
order of 1000°C in both air and vacuum. At the higher temperatures 
encountered in Case II, however, these insulations may experience 
shrinkage, decomposition and/or interaction with silicon-germanium 
alloys . 

7. If a protective gas is used inside the generator, the refractory foil 
offers the most advantageous insulation system. 

8. The graphite re-entry shell of the radioisotope fuel capsule in the 
MHW generator will operate at elevated temperatures. If the capsule is 

to operate in air it will be necessary to keep the shell from being oxidized. 


29 . 



In conclusion it may be stated that probable problem areas in the unsealed 
air operation of the MHW generator relate to the graphite re-entry shell of the 
fuel capsule and the foil insulation, if foil insulation that employs refractory 
metals is used in the generator. Both problem areas may be conveniently eliminated 
by sealing the generator and employing inert gases during air operation. If forced 
convection is used instead of the natural convection assumed in the present investi- 
gation, it may be possible to significantly decrease all generator operating 
temperatures. Although this procedure would somewhat alleviate some potential 
problem areas associated with the air operation of the "open" MHW generator, it 
would be by no means eliminate them . 
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TABLE I 


Characteristics of MHW-RTG Preliminary Reference Design 
As Described by General Electric 


RTG 

PERFORMANCE 

POWER 

VOLTAGE 
FUEL LOADING * 

t hot JUNCTION 

t cold JUNCTION 
CONVERSION EFFICIENCY 

GEOMETRY 
DIAMETER * 

LENGTH * 

CONFIGURATION * 

CONVERTER 

END CLOSURES 
LOAD STRUCTURES 
RADIATORS 

THERMOELECTRIC COUPLES 

NUMBER * 

MATERIAL * 

LEG DIMENSIONS * (Case I) 

LEG DIMENSIONS * (Case II) 

HOT SHOES 
INSULATION * 


148 WATTS BOM 

114 WATTS EOM (12 YR) 

30 VOLTS DC 

2000 WATTS (2-1-000 WATT HEAT SOURCES) 

1100°C (2012°F) BOM 

335°C BOM (635°F) 

7.4 % BOM 


11.85" ACROSS DIAGONALS (HEX) 
20 . 0 " 

6 SIDE STRUCTURE - NO FINS 


SPOKED TITANIUM 

6 TITANIUM LONGERONS 

BERYLLIUM -IRON TITAN. COATING 

e « 0.85 


288 (SERIES PARALLEL) 

SiGe (80% SILICON) 

N-LEG 0.8" LONG x .0188 DST 
P-LEG 0.8" LONG x .0208 DC 
N-LEG 1.2" LONG x .0238 IN^ 

P-LEG 1.2" LONG x .02 62 IN 
SiMo 

LINDE-TYPE FOIL, ZIRCAR, MIN-K 2020 


* Input data for performance calculations 
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TABLE 2 MHW-RTG Performance in Vacuum and Gas Environments 
for Thermoelement Length of 0.8 in. 


FOIL INSULATION 

GAS 

VAC 

HELIUM 

AIR 

ARGON 

KRYPTON 

XENON 

P_. Watt 
Max, 

138 

17 

69 

84 

106 

118 


137 


56 

80 

105 

117 

QE9I 

1065 

543 

810 

872 

939 

972 

T °C * 
X C, 

321 

282 

279 

278 

276 

275 

ZIRCAR INSULATION 

P., Watt 

Max, 

124 

24 

73 

88 

103 

111 

P 30V, Watt 

123 

- 

, 61 

85 

101 

109 

T„ °C * 

ri t 

1012 

591 

844 

885 

929 

950 

T °C * 

o. 

WSM 

282 


278 

277 

276 

MIN-K 2020 INSULATION 

P_. Watt 
Max, 

116 

35 

84 

94 

103 

108 

P 30V, Watt 

115 

- 

77 

91 

101 

106 

t h, ° c * 

986 

652 

860 

903 

929 

941 

T_ °C * 
O , 

323 

281 

279 

278 

277 

277 
■ ' 


* Temperatures are for maximum power points 
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TABLE 3 MHW-RTG Performance in Vacuum and Gas Environments for 

Thermoelement Length of 1.2 in . 


FOIL INSULATION 


GAS 

HI 

HELIUM 

AIR 

ARGON 

KRYPTON 

XENON 

P,, Watt 

Max, 

161 

28 

92 

114 

135 

145 

P 30V, Watt ■ 

155 

i 

91 

114 

132 

140 

__ o _ . 
T_ C * 

O. g 

1195 

641 

934 

1016 

1070 

1100 

T C * 

o. 

320 

282 

i 

277 

276 

274 

273 


ZIRCAR INSULATION 


P. _ Watt 

Max, 

135 

37 

101 

no 

| 123 

130 

P 30V, Watt 

134 

8 

101 

no 

122 

128 

T„ °C * 

xi 9 

1125 

700 

970 

1002 

1039 

1057 

„ o_ , 
T C * 

G, 

321 

281 

277 

276 

275 

275 


MIN-K 2020 INSULATION 


P,, Watt 

Max, 

123 

1 

51 

110 

105 

112 

n5 

P 30V, Watt 

122 

32 

110 

105 

112 

n4 

T„ °C * 

1102 

774 

1004 

987 

1006 

1015 

T ° c * 

O # 1 

1 

322 

280 

276 

276 

276 

276 


Temperatures are for maximum power points . 
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Fig . 2 Typical Generator Configuration 



GOLD SHOE 



POWER INPUT, WATTS 

Fig. 2. C-2 Min-K 2020 Insulated Converter Vacuum/Nitrogen Performance 
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Fig. 4 C-2 Min-K2020 Insulated Converter Vacuum/Nitrogen Performance 

, 37. 



POWER INPUT, WATTS 

Fig . 5 C~3 Foil Insulated Converter Vacuum/Nitrogen Performance 
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POWER OUTPUT, WATTS; EFFICIENCY 



POWER INPUT, WATTS 

Fig . 6 C-3 Foil Insulated Converter Vacuum/Nitrogen Performance 
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' >* ' * 



COLD SHOETEME. , C HOT SHOE TEMP*/ 



POWER INPUT, WATTS 

Fig. 7 C-3Z Zr 0_ Insulation Converter Vacuum/Nitrogen Performance 
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POWER OUTPUT / WATTS; EFFICIENCY % 



POWER INPUT, WATTS 


Fig. 8 C-3Z Zj-CL Insulation Converter Vacuum/Nitrogen Performance 
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EFFECTIVE THERMAL CONDUCTIVITY; WATT/cm °C x 10 



Fig . 9 Effective Thermal Conductivity in Vacuum and Nitrogen . 
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LOAD VOLTAGE, VOLTS HOT 



LOAD CURRENT, AMPERES 
Ficr-l 11 MHW Comparative Performance 44. 







LOAD VOLTAGE, VOLTS 



MHW - Generator Performance In Vacuum and Air 46 , 


Fig. 13 Foil'Insulated 


POWER OUTPUT, WATTS 


GOLD JUNCTION TEMP. , C HOT JUNCTION TEMP 



Fig . 14 ZrO^lnsuIated MHW - Generator Temperature In Vacuum and Air 
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THERMOELEMENT LENGTH = 0.8 in. 

T.. =311°K (100°F) 

AMB . 


Power Out 
In 

Vacuum 


Power Out 
In 
Air 


Load Voltage 

\ to — * 
vacuum 


Load Voltage, 
In -/ 
Air 


5 6 7 

LOAD CURRENT AMPERES 
MHW - Generator Performance in Vacuum and Air 


Fig, 15 Zr £L Insulated 


COLD JUNCTION TEMP. , G HOT JUNCTION TEMP 



CURRENT , AMPERES 


Fig. 16 MIN-K 2020 Insulated MHW Generator Temperatures in Vacuum and Air 
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'AGE, VOLTS 


Load Voltage 
In Vacuum 


THERMOELEMENT LENGTH = 0.8 In. 


t amb = 311 ° k (10 ° ° f) 


Power Output 
In Vacuum 


Power Output 
In Air 


Load Voltage 
In Air 


CURRENT, AMPERES 

Min-K 2020 Insulated MHW - Generator Performance in Vacuum and Air 


POWER OUTPUT, WATTS 





COLD JUNCTION TEMP. , C HOT JUNCTION TEMP. 



Fig. 18 Calculated MHW Generator Temp. - Foil Insulation 
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HOT JUNCTION TEMP. 
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THERMOELEMENT LENGTH = 1.2 in. 


T ... n = 311°K 
AMB 



LOAD CURRENT, AMPERES 

Fig. 21 CALCULATED MHW GENERATOR PERFORMANCE - Zr0 2 INSULATION 
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Fig .22 CALCULATED MHW GENERATOR TEMP. MIN - K 2020 INSULATION 
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Fig. 23 CALCULATED MHW GENERATOR PERFORMANCE MIN - K 2020 INSULATION 
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Fig .24 Thermal Conductivity of Gases 
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INSULATION CONDUCTIVITY, MW/cm 



GAS CONDUCTIVITY, MW/cm - °C 


Fig .25 Effect of Gas on Insulation Conductivity corresponding to average 
Temperature of the thermoelement. 
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(Air & Nitrogen) 







Cold Junction 


Peak Power 




TEMPERATURE 



Fig. 28 MHW-RTG Temperatures for Thermoelement Length =1.2 in. 
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APPENDIX A - AVERAGE THERMOPHYSICAL PROPERTIES 


The temperature dependence of various thermophysical properties re- 
quires the definition of average values for the properties to be used algebraically 
in thermoelectric generator performance equations . By assuming no radiation 
losses and no internal heat source/heat sink, the average thermal conductivity, 
Seebeck coefficient and electrical resistivity may be defined as follows: 

THERMAL CONDUCTIVITY 

The heat conduction equation may be written as 


q = constant = k(T) = <k>^^- 

dx i 


(Al) 


Integration of Eq. (Al) and the use of appropriate limits of integration yields 

T 

Hi j6 

<k> dx = <k> dx - <k>AT 


J k(T)dT =J <k> dx - <k> &Lj 


(A2) 


Upon rewriting, the average thermal conductivity may be defined as 

T 

r H 

<k> = srl k(I)dT 


(A3) 


SEEBECK COEFFICIENT 

The averaged Seebeck coefficient is obtained from its definition by 

integration 


H H H 

V = J S(T)dT =J <S>dT = <S>J dT = <S> At 


(A4) 


Upon rewriting, the average Seebeck Coefficient may be defined as 



ELECTRICAL RESISTIVITY 


The definition of average electrical resistivity requires the inclusion 
of a non-linear temperature profile. The average electrical resistivity is defined 
by 


<n> ~ = 




~~A P ft W ) dx 


(A6) 


Solution of Eq. (Al) for the differential dx, the substitution of this into Eq. (A6) 
and rearrangement yields 


< D > 



(A 7) 


The substitution of Eqs . (Al) and (A3) into Eq . (A7) finally gives 


H 


<p> = 


f p (T)k(T)dT 

2a 


I 


k(T)dT 


(Ad) 


Equations (A3), (A5), and (A8) serve to define the average thermophysical 
properties of the thermoelements . Providing the temperature limits of integration 
are known, the averaged properties are constant numeric values which may be used 
algebraically in the performance equations . 
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APPENDIX B - COMPUTER PROGRAM 


The computer program reproduced in this Appendix was used to perform 
all of the calculations reported in the present memorandum. The program is written 
in BASIC language and in its entirety is used to calculate generator performance in 
either vacuum or air. To calculate generator performance in vacuum operation, the 
convective coefficient, Ki, is set equal to zero. The effect of gases upon the thermal 
insulation is reflected in the value calculated for K4. The instruction GOSUB 720 
(line 690) is an output format sub-routine. The following definition of symbols 
applies to the computer program. 

Input 


Q1 

El 



A1 = A 

R 

A5 = A 

c 

K2 = <k > 
P 


K3 = <k > 
n 


A2 = a 

P 

A3 = a 

n 

N = N 


A4 = A 

e 

K4 = <kj> 

L = l 

R2 = p 

R3 = p 
*n 

3 =l<S n H+l<S p >| = 


total heat input, watts 

combined effective emissivity, fin effect and view 
factor 

radiative area, cm^ 
convective area, cm^ 

averaged p-type thermoelement thermal conductivity, 
watt/cm°K 

averaged n-type thermoelement thermal conductivity, 
watt/ cm 0 K 

p-type thermoelement area, cm 

2 

n-type thermoelement area, cm 
total number of couples 

2 

effective insulation area, cm 

effective insulation thermal conductivity, watt/cm°K 
thermocouple length, cm 

p-type thermoelement electrical resistivity, £J-cm 
n-type thermoelement electrical resistivity, Q-cm 
combined n-and p-type thermoelement Seebeck 
coefficient, v/°K 


R4 — 6 


electrical loss factor 



Output 


N1 = Ng 
N2 = n 

y -at c = 

K5 = K 


number of thermocouples in series 
number of parallel electrical branches 
cold -stack temperature drop 
convective heat transfer coefficient 


M = m 

T1=T C 
T2 = T 

V=V L 
P = P 

o 

e = v 

i = i 


cold junction temperature, °C 
hot junction temperature, °C 
load voltage, volts 
power output, watts 
efficiency, percent 
load current, amps 
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GET- GEN 

LIST 

GEN 

iO REM G1=HEAT INPUT* E 1 = £FF . EMISS I V ITY* A 1 = RAD I AT I V/E ARFA 
20 REM A5=C0NVECTIVE AREA* K2=T HFRM . C0ND 0F P-LEO 
30 'RFM K3=THERM. C0ND 0F N-LEC* A2=AREA 0F P-LEG 

40 REM A3= AREA 0F N-LFG* N=# 0F C0UPLES* A4= AREA 0F I NSUL AT I0N C EFF . ) 

50 RFM K A- INS. EFF. C0ND. * L=C0UPLE LENCT H* R2= P-LEG RES* 

60 RFM R3= N-LEC- RFS»* S=SEFBECK C0EFF * R4=ELECT . L0SS* N 1 =#SERI ES C0UP. 
70 RFM, N2 = # PARALLEL BRANCHES* Y= DELTA TC* K5=C0NVECTI VE C0EFF. 

H3 PRINT "ENTER VALUES F0R 0 1 * E 1 * Al * A5* K2* K3* A2* A3* N" 

90 INPUT G1*E1*A1* A5*K2*K3*A2*A3*N 

100 PRINT "ENTER VALUFS FOR A4*K4*L* R2* R3* S* R4*N1 *N2* Y*KS" 

HO INPUT A4*K4*L*R2*R3*S*R4*N1*N2*Y*K5 
120 F0R 1=1 T0 A . 

130 PRINT 
140 NEXT I 

'150 PRINT "PROGRAM 10 CALCULATE T/E C-ENFRAT0R PERFORMANCE" 

160 PRINT T AE (15)1 "HEAT INPUT="i01 
170 PRINT 
180 PRINT 

190 PRINT " M TC TH VL P0 EFF I” 

200 PRINT 
210 RFM 

220 FOR M=0 T0 2 STFP .2 
230 C0SUB 320 
240 NEXT M 

250 FOR M = 3 T0 10 STEP 1 
260 G0SUB 320 
270 NEXT M 

280 F0R M= 1 5 T0 60 STEP 5 

290 C0SUB 320 

300 NEXT K 

310 END 

320 02= . 95*01 

330 A=5.67F-12*E1*A1 

340 F3=A5*K5 

350 LFT c = 6 

360 X=0 

370 T = A» <-.25>*C.975*C 1 -E/1 005*01 ) t .25 
3P0 X = X+1 

390 F 1 =A*(T t 4-3 1 1 »4)+F3*CT-31 1 > t l .25- C 1 -E/ 1 GO ) * . 97 5*Q 1 
410 F2=4*A*T»3+ 1 .?5*F3* CT-31 1 ) t .25 

410 T=T-F!/F2 

420 IF ABS(F1/F2)< . 1 THEN 440 

430 GOTO 390 

440 RFM 

450 C = 7 + Y 

460 T 1 = C-273 

470 K= CK2*N*A2+K3*N*A3>* 1 /L+K 4* ( A4-N* A2-N* A3 ) * 1 /L. 

4P0 R 1= <R2/A2 + R3/A3)*L 

490 C=(N*St2)/C2*k4*Rl*( 1+M>r2) 

500 B=C2*M*C-K/GT*1/C 1 + P*M> 

51 0 D= CC2/ 0+(K/C)*C+CtF>*l/<l +2*M) 

520 H=CE+(Rt2+4*B)r .5)*l/2 

530 T2=H-273 

540 n=H-c 



550 V=CM/CH-M>)*N1*S*D1 

560 R5=<N1/N2)*R4*R1 

570 P= (M/ < 1 +M )tg)*(NI *S*D 1 ) 1 2* 1/R5 

580 F=P/Q 1 * ! 00 

590 IF X<4 THEN 380 

600 I = 0/(1+M>)*N1*S*DJ*I/R5 

610 ZU3=M 

620 ZC 2 3-T 1 

630 Z C 3 ] = T2 

640 ZL 4] = V 

650 Z C 5 3 =P 

660 ZC63 = E 

670 Z E 7 3 = I 

680 ZS="##.# #### #### ##.## ###.## ##.## 

690 C.0SUB 720 
700 RETURN 
710 STOP 

720 RFM ***** FORMAT ' ***** UTILITY SUBROUTINE ***** 
730 RFM ***** VFRSI0N 1 ***** 7/31/69 ***** 

7^0 RFM SUBROUT INF TO FORMAT A PROCRAM’S DATA 
750 LET Z2=Z3=Z4=Z5=Z7=Z8=Z9= 1 
760 DIM Y$C 1 0 3 # Z$C 723 
770 LFT Y $=”0 123456789" 

780 LET Z0=Z9-1 
790 LFT Z0=Z0+1 
POO IF Z0=LFN<ZS)+1 THEN 1310 
810 IF ZSCZO^ZOIs"#” THEN 880 
820 IF Z$CZ0*Z0+1 THEN 880 

830 IF Z£tZO.»ZO+l ] = ” + #” THEN 860 
840 PRINT ZSCZ0*Z03f 
850 GOTO 790 
860 LFT Z 4=0 
570 C0T0 790 
880 LFT Z=100 
890 LFT Z6=ZCZ23 
900 LFT Z 9=Z 0- 1 
910 LFT Z9=Z9+1 
920 IF Z SC Z9* Z9 3 = " • *' THEN 950 
930 IF ZS^:Z9,Z93 = ‘ , #• , THFN 910 
940 G0T0 990 
950 IF Z5#l THEN 990 
960 LFT Z 5=0 
970 LFT Z = Z9 
980 GOTO 910 
990 IF Z*100 THEN 1010 
1000 L FT Z=Z9 
1010 IF Z4=l THEN 1060 
1020 IF Z 6 > = a THFN 1050 
-"1030 PRINT "-"t 
1040 GPtTO 1060 
1050 PRINT ” "J 

1060 LFT Z6=ABSCZ6>+10»<Z-Z9-1> 

1070 ,,„F0 R Z1 = Z-Z0 TO Z+1-Z9 STEP -1 
1G60 IF ZS"CZ-Z1*Z-Z1 l#***” THEN 1130 
1090 PRINT 
1100 LET Z 3=0 



